A numerical study 

of wave-function and matrix-element statistics 
in the Anderson model of localization 

V. Uski 1 , B. Mehlig 2 , and R.A. Romer 1 

1 Institut fur Physik, Technische Universitat, D-09107 Chemnitz, Germany 
theoretical Physics, University of Oxford, 1 Keble Road, OX1 3NP, UK 

Received 6 October 1998, revised 16 October 1998, accepted in final form 23 October 1998 
by M. Schreiber 

Abstract. We have calculated wave functions and matrix elements of the dipole operator 
in the two- and three-dimensional Anderson model of localization and have studied their 
statistical properties in the limit of weak disorder. In particular, we have considered two 
cases. First, we have studied the fluctuations as an external Aharonov-Bohm flux is varied. 
Second, we have considered the influence of incipient localization. In both cases, the statisti- 
cal properties of the eigenfunctions are non-trivial, in that the joint probability distribution 
function of eigenvalues and eigenvectors does no longer factorize. We report on detailed 
comparisons with analytical results, obtained within the non-linear sigma model and/or the 
semiclassical approach. 

Keywords: Disorder, semiclassical approximation, wave function statistics 



1 Introduction 

Disordered quantum systems exhibit irregular fluctuations of eigenvalues, wave func- 
tions and matrix elements. The statistical properties of wave functions and matrix 
elements are of particular interest. Both are of direct experimental relevance. Fluc- 
tuations of wave-function amplitudes determine the fluctuations of the conductance 
through quantum dots. The effect of an external perturbation is described by matrix 
elements of the perturbing operator in the eigenstates of the system. 
In the metallic regime (which is characterized by a large conductance g ^> 1), wave- 
function and matrix-element fluctuations are described by random matrix theory 
(RMT) In Dyson's ensembles |] such fluctuations are particularly simple 

since the joint probability distribution function of eigenvector components and eigen- 
values factorizes j| and the statistical properties of the eigenvectors are determined 
by the invariance properties of the random matrix ensembles. One finds that the 
wave-function amplitudes are distributed according to the Porter-Thomas distribu- 
tion ||] . Non-diagonal matrix elements of an observable A are Gaussian distributed 
[gj around zero with variance tr^g. Diagonal matrix elements are also Gaussian dis- 
tributed, with variance a^ iag = ((3/2)a^ s where [3 — 1 in the Gaussian orthogonal 
ensemble (GOE) of random matrices and (3 = 2, 4 in the Gaussian unitary (GUE) 
and symplectic ensembles. The variance of non-diagonal matrix elements is essen- 
tially given by a time integral of a classical autocorrelation function Q and does not 
depend on the symmetry properties. 

Of particular interest are those cases where the fluctuations of eigenvalues and eigen- 
vectors are no longer independent. In the following we analyse two such situations. 
First, we consider the effect of an Aharonov-Bohm flux which breaks time-reversal in- 
variance and drives a transition from the GOE to the GUE. The statistical properties 



1 



of diagonal matrix elements in the transition regime between GOE and GUE were 
calculated in ]|, 0, [j| , those of level velocities in || 0, ^| . The statistical properties 
of wave functions in the transition regime were derived in ]l0| |. Here we compare 
these predictions with numerical results obtained for the Anderson model. Second, 
we study the influence of increased disorder on the statistical properties of wave func- 
tions. The question how the distribution of wave-function amplitudes deviates from 
the RMT predictions at smaller conductances g has recently been discussed in |pH , 
within the framework of the non-linear sigma model HHHI, and 

using a direct 

optimal fluctuation method |]l5| . Here we compare our numerical results for distribu- 
tions of wave-function amplitudes in the d = 3 dimensional Anderson model with the 
perturbative results of |p^| . 



2 The Anderson model of localization 

We consider the Anderson model of localization in d = 2 and d — 3 dimensions 
which is defined by the tight-binding Hamiltonian on a square or cubic lattice with 
N sites and unit lattice spacing 

H = \n)e n (n\ + \n)t nm (m\ , (1) 

and |n) represent the Wannier state at site n. The on-site potential e n is assumed 
to be random and taken to be uniformly distributed between —W/2 and W/2. The 
hopping parameters t nm connect only nearest-neighbour sites (and t = 1). In the 
presence of an Aharonov Bohm-flux ip, the hopping parameters acquire an additional 
phase. The flux tp is measured in units of the flux quantum tpo = he/ e and we define 
4> = f/fa- The presence of an Aharonov-Bohm flux (f> breaks time-reversal invariance. 
We have determined the eigenvalues E a and the eigenfunctions ?p a (r) using a modified 
Lanczos algorith m The statistical properties of eigenvalues, eigenfunctions and 
matrix elements A a p = (il; a \A\ipp) of an observable A depend on the strength 
of the disorder potential. In the metallic regime (which is characterized by a large 
conductance g) these fluctuations are described by RMT on energy scales smaller than 
the Thouless energy Ed = gA, where A is the mean level spacing. This is the regime 
we consider in Sec. 3, while Sec. 4 deals with g" 1 corrections to the distribution of 
wave- function amplitudes. 

3 The GOE to GUE transition 

In this section, we discuss numerical results for the fluctuations of matrix elements, 
level velocities and wave functions in the transition region between GOE and GUE. 
We have calculated the smoothed variances 

C v (e,0) = {\d v (E^-t)\ 2 ) E C m {e,<p) = {\d m {E,fre)\ 2 ) E (2) 

where d v , m (E, <p\ e) are the fluctuating parts of the densities 

d^E^-e) = Y-gf^E-E^)} (3) 

a 

d m {E,cp-€) = J2 A ^ S ^ E - E «(W ( 4 ) 

a 
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Fig. 1: Velocities and matrix elements in the transition regime between GOE and GUE 
for the 27 x 27 Anderson model with W = 1.7 (symbols), compared to the corresponding 
semiclassical expressions (lines). The numerical values of e are 0.158 (v)> 0.224 (x), 0.316 
(<), 0.447 (o), 0.631 (□), 0.891 (O), 1.26 (A), 1.78 (+), 2.51 (*). 



with 6 e (E) = (V27re)" 1 exp(- J B 2 /2e 2 ). It is assumed that (A aa ) = (dE a /d(j>} = 0. 
C m was calculated exactly in ||. Here we compare with semiclassical expressions 
for C m |J f?j and C v Q obtained within the diagonal approximation and valid for 
A < e < gA (see also ||). In the case of C m we considered the matrix elements 
of the dipole operator A = x. The numerical data were obtained by averaging over 
69 realisations of disorder in the 27 x 27 Anderson model with W — 1.7 and over all 
states in the energy interval [—3.4,-1.9]. The two-dimensional case was considered 
in order to be able to obtain a good numerical accuracy at each flux value in a 
tolerable computing time. In general we observe good agreement. For a more detailed 
discussion of the results see jj| . 

In Fig. 2(a) we show corresponding results for the distribution function of wave- 
function amplitudes which is defined as follows 

f(t) = A(^25(t-\ip a (r)\ 2 N)S(E-E a ^J . (5) 

Within RMT one obtains for f(t) in the limiting cases of GOE and GUE 

/ P G T OE (t) = _Lexp(-f/2), (6) 
/ P G T UE (t) = exp(-t) (7) 

(Porter-Thomas distributions |3|]). An expression in the transition region between 
GOE and GUE was derived in Jl0|. We have computed f(t) numerically for several 
values of <j) in the 13 x 13 x 13 Anderson model at W — 0.5 (in the metallic regime). 
Our results are shown in Fig. 2(a), where we have plotted 2y/tf(t) versus y/i. The 
results are in good agreement with the analytical formulae (6), (7) and those given in 

& 

It was predicted in |6| that the distribution of velocities ceases to be Gaussian in the 
transition regime between GOE and GUE. The deviations are small, however, and 
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Fig. 2: Wave-function statistics in the 13 x 13 x 13 Anderson model, (a) GOE to GUE tran- 
sition of the distribution function in the metallic regime. The Porter-Thomas distributions 
are shown as dashed lines. The predictions of [[lo| are show as solid lines, (b) Deviations 
Af(t) from the Porter-Thomas distribution for the orthogonal case, and several values of 
disorder strength, W = 1,2,3,4 and 5. The results are averaged over the energy interval 
[—1.7,-1.4]. The dashed vertical lines denote the zeros of the first order correction term in 
Eq. The solid lines are fits according to Eq. (^). 



we have not been able to reduce the statistical errors of our numerical results to an 
extent that a meaningful comparison with the results of || becomes possible. 

4 Deviations from RMT 

Within the non-linear sigma model it is possible to derive g _1 -corrections to RMT. 
For distributions of wave- function amplitudes this was done in fll3| , where in the 
orthogonal case one obtains 

/(*) = /pt° E («) [1 + ^g- 1 (3/2 - 3t + t 2 /2) + O {g- 2 )] . (8) 

In the case discussed in Ref. |L3|], ~ LI I where L is the linear dimension of the 
system and I is the mean free path. Eq. (j3|) is valid provided t <C sj g/cid- 
We have performed numerical simulations in the d — 3 Anderson model, using dif- 
ferent values of disorder, and have computed the distribution f(t) by averaging over 
400 realizations of disorder. The deviations from the Porter-Thomas distribution, 
= f{t)/ fpT E (t) ~ 1' 1S shown in Fig. ||(b). In all cases, the deviations exhibit a 
characteristic form: the probability of finding small and large amplitudes is enhanced, 
while the distribution function is reduced near its maximum. We find that A/(i) can 
be fitted using Eq. (g). At increasingly large disorder, deviations occur at lower values 
of t. In all cases, however, the zeroes of Eq. (||) are well reproduced. Since at weak 
disorder I ~ W~ 2 one might expect a^/g ~ W A . In the present case, however, we 
find a-ijg ~ W 2 . In order to resolve this discrepancy, more accurate numerical data 
at small values of W are needed. Corresponding results have been obtained for the 
unitary case. 
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5 Conclusions 



We have analyzed RMT fluctuations of matrix elements and level velocities in the 
d = 2 dimensional Anderson model, in the transition regime between GOE and GUE 
and have found good agreement between the predictions of RMT and our numerical 
results. For the distribution of wave-function amplitudes in d = 3 we have studied 
deviations from RMT, in the form of g - ^corrections, as suggested in [pi. Our nu- 



merical results can be fitted by the expressions derived in 13 , the dependence of 



the fit parameter on W however, differs from what might be expected. In this 
context it will be very interesting to determine the tails of the distribution functions 
and compare with the predictions in p2j, ftM and [15| . 
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